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An improved vortex filament-vortex core method for predicting aerodynamic 
characteristics of slender •''lags with edge vortex separation has been deve- 
loped. Semi-empirical but simple methods are used to determine the initial 
positions of the free sheet and vortex core. Comparison with available 
data Indicates that: (1) the present method is generally accurate in pre- 

dicting the lift and Induced drag coefficients but the predicted pitching 
moment is too positive; (2) the spanwlse lifting pressure distributions 
estimated by the one vortex core solution of the present method are 
significantly better than the results of Mehrotra’s method relative to the 
pressure peak values for the flat delta; (3) the two vortex-core system 
applied to the double delta and strake-wing produce overall aerodynamic 
characteristics which have good agreement with data except for the pitching 
moment; and (4) the computer time for the present method is about two thirds 
of that of Mehrotra's method. 
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2. Introduction 


It la w*. .1 known that at modarata to high anglaa of attack, tha flow 
aaparataa going around tha laadlng-adga of a aharp-adged wing. Tha 
aaparatlon producaa rollad-up vortax ahacta abova tha wing aurfaca which 
aaanata from tha local laadlng-adga. Tha concantratad vorticlty abc.va 
tha wing lnducaa high auction praaaura over a larga portion of tha uppar 
aurfaca of tha wing, raaultlng In a nonllnaar contribution to tha wing 
loading. In tha analyala and daalgn of alandar high apaad airplanaa 
flying at off daalgn conditions, a datailad knowledge of tha particular load- 
ing on tha wing aurfaca la raqulrad to aatlmata aarodynamlc performance 
and atructural loada. 

Thara hava baan several dlffarant theoretical approaches In aolving 
tha vortax flow aaparatlon problam ainca tha aarly 1950 'a. In general, 
they can ba claaalfiad aa conical and non-conlcal flow mathoda. 

In the flrat category, baaed on the Blander wing approximation, the 
aolutlons are obtained from a simplified governing flow aquation by 
neglecting the variation of longitudinal velocity. Therefore, the problam 
la solved In the so-called cross-flow plane. Legendre (ref. 1) assumed 
that the rollad-up vortax sheets could be replaced by a pair of concentrated 
cores over the wing and Inboard of the leading-edge. Two boundary 
conditions were Imposed such that velocities along the leading-edges 
must be finite and the total force on each Isolated vortex must be zero. 

Later, Brown and Michael (ref. 2) used feeding cuts to connect the line 
vortices and the wing leading-edges. Again, the boundary conditions require 
tangency of flow on the wing surface, finite velocity at the leading-edges, 
and zero resultant force on the vortex-plus-cut combination. Mangier and 
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and Salth (ref. 3), aodlfylng the nodal In rafaranca 2, aaaunad that tha 
vortax ahaaca ara ahad tangentially at leading-edges and roll up an aplral 
vortlcaa. Tha boundary condltlona of caro pressure discontinuity acroaa 
tha aaparatad vortax ahaat and zero total forca on tha vortax shaat and 
faad cut arc laposad. Finally, Smith (raf. 4) davalopad tha fardel of a 
sagnantad feadlng vortax shaat by improving tha numarlcal procadura 
in rafaranca 3. Tha main ahortcomlngs of all thesa modala daacribad 
above ara llmltad to conical flow. Tha aaaumptlon of slandamass loads 
to a thaory which is lndapandant of Mach numbar. Thus, tha thaory 
braaks down at low Mach numbars in tha tralllng-adga regions , bacause 
tha traiilng-edge Kutta condition is not satisfied. It braaks down also 
whan tha aspect ratio is greater than ona. 

In the second category, tha conical flow assumption is removed so 
that the Kutta conditions ara satisfied at both leadlng-and t railing-edges. 

In the earlier stage of development, tha classical lifting surface theory 
was extended by Including a relatively simple representation of the leading- 
edge separation. Bollay (ref. 5) suggested a method for low aspect ratio 
rectangular wings in which all the vortlclty is shed from the wing aide-edge 
and is assumed to be along straight lines making an angle of a/2 with respect 
to the wing planform. Gersten (ref. 6) and Garner and Lehrlan (ref. 7) 
extended Bollay' s model to cover wings of arbitrary shape. Vortlclty la 
assumed to be shed from the leading-edges at an angle of a/2 to the wing 
plane. The wing was replaced by a series of lifting elements. A Multhopp-type 
lifting surface method (ref. 8) was used to solve the integral equations in 
reference 7. 

Nangia and Hancock (ref. 9) combined the conventional lifting surface 
theory with the detached flow model of Brown and Michael (ref. 2). 
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In the latter modal, two concentrated coraa 11a above the wing, connected to 
the leading-edges by feeding cute. The trailing wake conelete of vortlcee 
ahed from the tralllng-edge and the continuing concentrated vortlcee. To 
determine the separated vortex strengths, their position and the wing vor- 
tlclty distributions, the flow tangency condition on the wing surface and the 
trailing-edge Kutta condition are applied at selected collocation 
points, together with the conditions of zero loading on the vortex 
and cut. The leading-edge Ki ta condition it automatically satisfied 
by the choice of the loading function. 

When the wing loading distribution is not required, Polhamus' method of 
leading-edge suction analogy (ref. 10) has been found to be simple and accurate. 
He assumed that the overall lift force can be estimated from the sum of 
potential and vortex lift contributions. The potential lift Is defined as 
the lift due to potencial flow about the wing with zero leading-edge suction. 
Vortex lift is caused by large suction forces induced by vortices 
in the separated flow, acting on the upper surface of the wing, and 
is assumed to be equal to, in magnitude, the leading-edge suction force in 
the attached potential flow. 

The detailed literature survey of the above flow concepts were 
reviewed by Matoi (ref. 11) and Parker (ref. 12). 

Since 1970, several methods have been developed by including a more 
precise representation of the leading-edge separation. Two popular methods 
under the second category are of the doublet-panel type (refs. 13, 14 and 
15) and the free, vortices (refs. 16, 17 and 18). In the method of doublet- 
panel type, such as the free-vortex-sheet method of Johnson, et al. (refs. 

13 and 14), the wing surface and free sheet are replaced with doublet 
panels that have biquadratic strength. The free sheet, whose shape and 
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position must be determined by the iteration, ie attached ur. >er all conditions 
to a fed sheet. The fed sheet approximately represents the vortex core 
affect and moves with the free sheet. The entire set of doublet utrengths 
is determined simultaneously during the Iteration process. Hoeljmalters 
(ref. 15) extended the method of Johnson, et al. (ref. 14) by including the 
entrainment effect of rotational cores. In this model, the spiralling free 
shear layer springing from the leading-edge is replaced by a free vortex sheet 
which is terminated by a combination of a feeding sheet and a discrete line 
vortex and sink. The strength of the sink is related to the entrainment of 
the rotational core and derived from seml-empirlcal arguments. The combi- 
nation of vortex/slnk and feeding sheet sustains no resultant force and pro- 
duces the condition needed for determining the vortex core position in the 
cross-flow plane. From references 13 to 15, the panel method has been shown 
to predict accurate results for flat delta wings. However, 
the computing time involved is lengthy and the methods are generally re- 
stricted to simple geometries. 

In the method of free vortices, Mook and Maddox (ref. 16) have used an 
extension of the vortex lattice method to account for leading-edge separation. 
In addition to standard vortex lattice, a series of kinked vortex lines are 
assumed <:o be shed from the leading-edge. An iterative procedure is used 
to determine the strength and local flow directions with appropriate boundary 
condition. In the Kandll method (ref. 17), the flow is modelled in a 
manner similar to that of Mook and Maddox (ref. 16), but included wing-tip 
separation. In this model, the bound elements of the vortex lattices were 
unswept, and the wake behind the tralling-edge and tip vortex elements 
were force free. The qua^i-vort ex-lattice potential flow method of Lan 
(ref. 19) has been extended by Mehrotra (ref. 18) to include vortex flow 
effects. This is done by modelling with discrete trailing vortex filaments 
emanating from the leading-edges. The leading-edge Kutta condition is exactly 
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satisfied snd partial span vortex separation Is allowed. The results of Meh- 
rotrs's model (ref. 18) Indlcste that the prsdlstsd spanwlse pressure distri- 
bution Is more diffused than the experimental data show. This Is because the 
effect of a concentrated vortex core In the real flow Is not well represented 
by a number of free vortices. From the published results, It Is clear that 
none of the existing free-vortex methods Is cspable of predicting accurstely 
the lifting pressure distribution over the wing surface within a reasonable 
amount of computer time. 

The main objectives of this investigation are : (1) to reduce the computer 

time through the use of a pair of concentrated vortex cores above the wing sur- 
face, which should also improve the pressure prediction of Mehrotra's model 
(ref. 18); (2) to allow the effect of side-edge vortex separation; and (3) to 
extend the method of reference 18 to more realistic geometries. Configurations, 
which are applicable to real airplane geometries Include features such as 
camber, strokes and leading-edge vortex flaps. All these configurations will 
be treated in this report. 

In the following, the aforementioned improvement will be described. First, 
the leading-edge vortex system is superimposed on the regular quasl-vortex- 
lattlce grid by using segmented vortex filaments, as shown In reference 18. 
Second, the side-edge vortex system Is obtained by extending all the bound 
elements to the tip chord In the last vortex strip. Third, the leading-edge 
diffused vortex elements. In reference 18, are allowed to merge into a con- 
centrated core after two iterations. 
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3. Description of Theoretical Formulation 


Mehrotra's method (ref. 18) is e diffused vortex-filament model of the 
vortex core but serves as a basis for the present method, in th< present 
analysis,' the basic assumptions are: (1) tho wing is represented by a bound 

vortex sheet across which there exists a pressure difference; (2) the sep- 
arated flow along l<Mdlng~edges is represented by force free leading-edge 
vortex elements which feed vortlclty into a concentrated core. In the present 
method, the Quasl-Vortex-Lattlce method (QVLM) (ref. 19) is used to simplify 
the induced velocity expressions due to the wing bound vortex sheet. 

During the iteration process, the following boundary conditions are 
imposed: (1) the wing surface must be impermeable; (2) Kutta conditions are 

imposed along the leading- and trailing-edges of the wing; (3) in the diffused 
vortex- filament model, the force free condition is applied on the leading-edge 
vortex filaments and trailing wa^e elements; and (4) in the core model, the 
force free condition is applied on the free vortex elements, concentrated 
core and trailing wake elements. 

This is a non-linear problem because the strength of the wing bound vor- 
tices and free element vortices, the locations of the leading-edge free 
vortex elements, and the concentrated core are unknown. 

3.1 Mathematical Model of the Wing 

The wing lies in the x-y plane and the x-axis is taken to be streamwise 
along the wing center line. The origin of the right-handed rectangular coor - 
dinate system is assumed to be at the moment reference point of the wing. 

3.1.1 Wing Bound Element and Trailing Vortices 

According to the QVLM, the semi-circle method is used to determine the 
control station locations and the vortex strips in the spanwlse direction, and 
the wing control points and bound elements in the chordwlse direction. The 
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locations of bound- and t railing-vortex elements ara aa ahown in Fig. 1. The 
x-locatlon of bound element a la given by 


*i * x t + 7 <1 - co8 < "" 2 ' ^"" ^) ) ♦ 1 " 1 * N 

where la the leading-edge x-coordlnate, c ia the local chord and N is 
the number of bound elements in the chordwlae direction. The spanwlse lo- 
cations of trailing vortices are given by 


yj - j (1 - cos ( 21 - jl - M ), j - 1, M (2) 

where b is the wing span and M is the number of legs of trailing vorticity, 
which is one higher than the number of spanwlse vortex strips. The locations 


of control points are given by 

x__ - x„ + ^ (1 - cos (jfir)), k - 0, 1, N 


CPt 


j 


cp 


j 


^ (1 - cos (^ir)), j - 1, (M - 1) 


(3) 

(4) 


where x^ - and Cj are the leading-edge x-coordinate and local chord at y^ , 
respectively. 

For the st rake-wing configuration, there are two sections in the span- 
wise direction. Each spanwlse section is divided into vortex strips by the 
semi-circle method. The edges of the vortex strips in each section are ob- 
tained through the following relation 

b. 




'k-i + T (1 ' cos \ 


and y control points are given by 
. b 


CP. 


y k-i + T (1 ‘ cos ( ^ ir)) ’ 1 “ l ' — V 1 


th 


(5) 


( 6 ) 


where b^ is the width of the k spanwlse section and is the number of 
spanwlse strips plus one in section (see Fig. 2). 
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y, , arc defined as 
7 k-l 





3.2 Model of Free Vortex Element System 

3.ZJ Leading-Edge Vortex System Geometry 

The le8ding-edge vortex system is superimposed on the regular quasi- 
vortex-lattice grid. A typical leading-edge vortex element is shown by 
points A through J in Fig. 3. These points are connected by a series of 
short straight segments. A typical set of initial locations for these 
segments is indicated with dashed lines and the locations after two Iterations 
are given by solid lines. 

The segments of the vortex element have the following characteristics: 

(1) Initially, points A through E lie along a wing trailing vortex 
filament with point A being one root chord away from the 
trailing-edge (point D) in the downstream direction, and the line 
segments between A and D are parallel to the axis of symmetry. The 
segments B-C and C-D are 0.1 c^ long. The length of each line segment 
between A and B remains unchanged as the solution progresses. Segment 
B-C is allowed to move only in the vertical direction, whereas segment 
C-D is fixed in the wing plane because of the trailing-edge Kutta 
condition. Segment D- E is also fixed in the wing plane. 

(2) Points E, F, G and H lie in the wing plane. The location of segment E-F 
is £ head of the wing first bound element and is given by 






X 

E 


( 7 ) 


a f , , tt % v 

f ' 1 • coa < ' f (^1) ' ' 


%v ' \ + -r (1 " co ‘ < 7?fey> > ( 8) 

The eegDente F-G and G-H are of the same length end point G lies on the leading- 
edge. These segments are fixed In the wing plane to satisfy the leading-edge Kutta 
condition. 

(3) The initial location of point T Is a function of a and Is given by 
(ref. 18) 

" X F (9) 

■ v_ (10) 


0.1 Cp Tan (22.5 - 0.5 a), for 


15 


( 11 ) 


or 


Z I ■ 0.1 Tan a, for a > 15° (12) 

Where is the root chord and a Is the angle of attack. Initially point 
J is one root chord away from the traillng-edge. The line connecting 
I and J may be divided into segments of equal length which lie in a 
plane parallel to thex-z plane. The initial height of these segments is 
approximately 0.1 c^ above the wing plane. In the final converged location, 
all segments between point H and J (free vortex elements, concentrated core 
and trailing wake elements) are aligned in the local velocity direction at 
their mid-points to satisfy the force free condition. 

(4) The semi-infinite segments from point A to infinity and J to infinity are 
straight and are parallel to the undisturbed free stream direction. Since 
the far wake of the leading-edge vortex element has very little effect 
on the wing aerodynamic characteristics, the far wake is assumed to be semi- 
infinite and parallel to the free stream. 
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3.2.2 Side Edge Vortex System Geometry 


To establish the side-edge vortex system as shown In Fig. 4, the leading- 
edge vortex element and wing bound element of the last vortex strip Is modified 
ss follows: 

(1) Extend the leading-edge bound element to the side-edge (point D). Points 
D, E, F, G and H lie In the wing plane. The segments F-G, G-H, H-l end 
I through J are given in the seme manner as in the preceding section. The 
extended leading-edge bound element consists of segments D-E and E-F. 


The location of each point Is given by 
“Til 


*D 


Tip 


+ -y-E ( 1 - Cos ( $-+!)> > 


c. 


X F " X i F + 1 (1 " C ° B ( 2 (N+1)* } 


X E " 2 + X F^ 

(2) The initial location of each vortex segment along the trailing leg of 
the extended leading-edge bound element, point A through D, is given by 


(13) 

(14) 

(15) 


x J+1 - Xj + Ay cos j 


(16) 


^j+1 y Tip 


Z J+1 - z j + ^ 2 


(17) 

(18) 


where 

Ay ■ 0.93 x (width of the last vortex strip). 

(3) Extend the regular wing bound element to the side-edge. (Points P, Q, R, 
S). Points P, Q, R and S lie in the wing plane and at the side-edge. The 


location of each point is given by 

. ^Tip ,, ,21 - 1 xv , 

x - x £ + — y d ” cos ( — 2N — 1 " N 

Tip 


y i ’ y Tip 


(19) 

( 20 ) 
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where x. la che leading-edge x coordinate ac Che aide-edge of Che wing; 

Tip 

che aide-edge local chord length. 

(4) Initially, all Che segments between polnta D and A, P and T, Q and V, etc. 
lie In a plane parallel to x-z plane at the wing aide-edge. Pointa A, T, 
U, V and W are one root chord away from the tralllng-edge. The initial 
aide-edge vortex aegment along polnta P through T (or Q through U, etc.) 
ia given by, 


x 1+1 ■ x t + Ax coe j 


( 21 ) 


i+1 


' 1+1 


'Tip 


z^ + Ax sin j 


( 22 ) 

(23) 


where Ax la assumed to be half length of the side-edge chord. 


(5) The semi-infinite segments from points A, J, T, U, etc. to infinity are 
straight and are parallel to the undisturbed free stream. 


3.3 Induced Velocity Computet ion 
3.3.1 Induced Velocity Due to Wing 

In the quasl-vortex-lattlce method (ref. 19), to satisfy wing boundary con- 
dition, the continuoua vortex dlatrlbution over the wing ia replaced by a quasi- 
contlnuoua one, being continuoua chordwiee but stepwise conatant in the apanwlae 
direction. Thue, the wing surface can be divided into a number of vortex 
stripe with the associated trailing vortices (Fig. 1). In any strip, consider 
a vortex element y dx with an arbitrary direction t (Fig. 5). The induced 
velocity due to all bound elements in the i C ** strip is given by (ref. 19) 


5. <«') \ ,2 ( 

ii 4lt Xj It I .. J I I l 


V 


k l n | a * x t' | 2 |fc'| I a ' | 

and due to the associated trailing vortices by 

(t) - 8 2 f*t y (*') (/» (*1 - *) x d t . , ) 
u TUx - • 1 


ai __J .T'dx’ (24) 


'6 


(25) 


The transformation, x' ■ x^ + c(y) (l-cos9 )/2, reduce equations (24) and 


(25) to 


and 


$ ( $) . g 2 c(y) / \ (9) y (9) Sin 0 d0 

ij 8if o 


* ,ts . g 2 c(y) f* 2, ( 9) Y ( 9) Sine d 0 
V i- W 5 TT *0 


(26) 


(27) 


where 


a x I / 

1 2 ' 1 -‘” 


- a 


!!• *H 2 ‘ |«’| |S’| 


-)• t 1 


« 2 w) -jt ^ 


and c(y) s ^ - x a, etc. are defined in the List of Symbols. 
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The total lnducad valoclty dua to tha i strip of vortex distribution 


is givan by 


ft t (ft) • J 0 ft x (e) y (e) Sine c 

2 

+ 8^ y > f *2 (8) Y <0) Sln ° dfl 

0 

2 Tt 

' ^2 (e) Y(fl) Sir,e de 


Whare tha first term is dua to bound element, sacond dua to laft lag of trailing 
alement and third dua to tha right lag of trailing alamanta. Tha abova lntagrals 
are reduced to finite suns through the mid-point trapezoidal rule (raf. 19): 


ft t (ft) - ---gl 7 — r (ft, + ft; - ft;’ ) Y k Sin 0. 

8N k-1 h 2 k Z k k k 

where 9. » (2k - 1) n/2N and locations of bound elements are given by 


x. the leading-edge x -coordinate at y, (left leg) 
1“ 1 

x. ■ the leading-edge x-cordinate at y, (right leg) 


c^ ■ chord length at 
C2 ■ chord length at y^ . 

The location of control points and trailing vortices are described in equations 
(1) through (6). 

Thus, the Induced volocity due to all vortex strips of the wing can be 


written as 


M-l 

ft (ft) - E ft, (ft) 

* 
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3.3.2 Induced Velocity Due to Free Vortex Eleaent Syetea 


The leadlng-edg • free vortex eyetea cone late of M - 1 distinct vortex 


elements. Each eleaent nay have different number of snail vortex segments, 
th 

Assume that the 1 set has L small segments. The Induced velocity at a given 
point (x,y,z) due to J C ** segment of i** 1 element Is given by (ref. 19) 


<*> 


b 2 i\ 


a x 




):• * tf \ p \ |;*| 


-) .t« 


(34) 


where 


I ■ x t + y J + z ft 
a ■ (Xj-x) t + (y^ - y) ? + (z^ - z) 1c 
5 ■ (x j+1 - x) t + (y J+1 - y) T + (*j +1 - *) ft 
ft ■ < x j + i _ Xj) t + (yj +1 - y) T + Czj +1 - Zj) ft 

a' * (Xj - x) ft + $ (y^ - y) ft + S (z^ - z) ft 

ft' " (x j+l “ * + 6 (y j+l~ y > ft + 6 (*j +1 - z > & 

ft' - (Xj +l - Xj) ft + 0 (y J+1 - yj)J + 8 (*j +1 - z j> ^ 

T* ■ leading-edge bound vortex strength of i** 1 strip. 

The subscripts J and j + 1 correspond to the end point of j 6 * 1 segment. The 
Induced velocity due to l t ^ element can be written as 


ft* (ft) •'ZZ ft ** (ft) (35) 

1 J-l 13 

Therefore, the induced velocity due to the entire leading-edge free vortex 
elements Is 


ft (ft) 


M-l 

z 

1-1 


u. 


(ft) 


(36) 


3.4 Formulations of Wing Boundary Condition 

The bound elements and the corresponding control points of the wing 
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surface ars numbered from ths leading-edge to the tralllng-edge and from 


tha root to tha aide-edge. The flow tangency condition can ba written aa 


KjI'vK) (r v ■ 


Coaa (~| )* - Sina 


a lk + Coaa (~-) k - Sin a 


4 - 

N. 


( 3 ?) 


whara N is tho total nuabar of wing bound elements; N. tha nuobar of 

leading-edge vortax alananta; tha lnducad downwash at tha l** 1 control print 

cf wing due to a unit horsahoe vortax density at j; B^ k tha induced downwash at 
th th 

the 1 control point due to the k leading-edge vortax element of unit 

th 

strength; Yj the vortex density of tha j bound element; tha strength of 
the leading-edge vortax element; (dz/dx)^ tha camber slope at the l** 1 
control point and a the angle of attack. According to equations (52) and (53) 
of reference 19, the leading-edge thrust coefficient Is related to the normal- 
wash on the leading-edge by the following equation 


a^ k ■ (Induced normalwaab - Cos a (^“) k + Sin <Ol at ^ e 


(38) 


which leads to the last expression on the right hand side of equation (37). In 


the expression, a £k Is defined as (ref. 19) 


Ik 


N J Tan 2 A + 0 2 ( pL 


2 c„ Coo A 


where 


ir Jl - M 1 Cos^ A 

” 00 




N ■ number of chordwlse bounded vortex elements 
A ■ leading-edge sweep angle 

M ■ free stream Mach number 

00 



, • sectional leading-edge thrust coefficient of the k fc ^ strip 

C k 

For complete leading-edge separation cases, the sectional leading-edge suction 
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1> tiro and ao is a 


ik* 

3.5 Forca Fraa Condition 

Tha vortax sagmants of each fraa alaoant abova tha wing aurfaca and Cha 
trailing waka ara to ba allgnad In tha dlractlon of local velocity vector cal- 
culated at thalr mid-point a. Conaider tha i C ** aegment of a vortax element. 

Tha coordlnatea of lta and polnta ara given by (x^, y^, z^) and (x t+1 , y^> 
z^). Aaauma that tha Induced velocity at tha mid-point of thia aegment at 
a given iteration la given by (ut + + «&). In the following sections, tha 

formulatlona of force free condition will ba derived in tha dlffuaad vortax 
filament and core modela reapectively . 


3.5.1 Formulation of Forca Free Condition of Free Vortex Elemanta in The 
Dlffuaad Vortex Filament Model 

For the i 6 * 1 aegment, the new locatiou of the (i+l) C ^ and point can be ob- 
tained by aligning the aegment in the direction of local velocity vector. 

That is 


x i+l 

*i+l 

z i+l 


y i 

z i 


♦5 


(40) 

(41) 

(42) 


where 

U - (u 2 + v 2 + w 2 ) 1 * 

Aa - ((x i+1 - x t ) 2 + (y 1+1 - y t ) + (z i+1 - z ± ) ] ** 

Before equations (40) - (42) are used, the following points should be considered: 

(1) The length of each aegment is to be preserved. 

(2) The free vortex segments above the wing should not come too close to 
the wing surface to avoid numerical difficulty in the present inviscid 
theory. 

(3) The adjustment of the location of each aegment to satisfy the force 
free condition should be such that it does not result in numerical 
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fluctuation*. 


Bated on the above cone ide rat Iona, equation* (40) to (42) will be modified aa 
follows : 

Consider the same 1 th segment. If this segment moves "a" percent only 
according to the velocity computed at its mid-point, the equations (40) to (42) 
can be modified to be 

(43) 

(44) 

(45) 


Ay' - a As + (1 - a)(y i+1 - y t ) 

Az' - a As + (1 - a)(* 1+1 - s^) 

Ax' - [as 2 - (Ay') 2 - (Az 1 ) 2 ] ** 


It follows that 


i+1 


i+1 


'i+1 


x A + Ax' 
y t + Ay' 
z^ + Az' 


(46) 

(47) 

(48) 


The free elements of the leading-edge vortex system have been restricted not to 
come any closer than a minimum specified height to the wing surface, which is 
given empirically by Mehrotra (ref. 18) as follows 


min 


» 0.1 c_ Tan (22.5 - 0.5a) , for a <_ 15 


z . • 0.1 c., Tan a 

min R 


, for a > 15 


(49) 

(50) 


If z. ,, is less than z . , it is then set z... equal to z . and Az' is re- 
t+l min l+i min 

calculated by using 


Az' ■ z . - z. 

min i 


(51) 


This restriction is needed because whenever the free elements are close to the 
wing surface, they induce large velocities on the wing and vice versa, which 
makes the locations of free elements fluctuate. In the real flow, at small 
angles of attack, the leading-edge vortex system is weak and diffused. The 
effect of the free vortices is artificially reduced by increasing z fflin as the 
angle of attack is decreased below 15 degrees. Note that equations (43) through 
(48) described above are mainly used to (1) establish the leading-edge vortex 
filament system in the first two iterations; and (2) find the new position of 
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--- v--M. j j g;„ | f , r - |n . . . 


trailing wake elements and aide-edge vortex ayatam in tha antira itaration 
procaaa. 


3.5.2 Foroulation of Forca Fraa Condition of Frea Shaat Vortax Elamant in Tha 


Cora Modal 

For tha i C ^ aagmant along a fraa vortex elamant between tha leading-edge 
and concentrated core, the initial alopes are determined by the total in- 
duced velocltiea at the mid-point of thia aegment 

(52) 


<£*)' 

v Ax ; 


( M)’ 

'Ax' 


v 

u 

w 

u 


.th 


The new slopes of the i segment of a free vortex element are given by: 
- d* (£) + (1 - d*)(&)' 


<&>* - d * <s> + (1 - 


(53) 

(54) 

(55) 


Similar to the parameter "a" in equations (43) and (44) of section 3.5.1, the 
relaxation parameter of free sheet, d, is used in equations (54) and (55) in 
the core model, to avoid any rapid variation of free sheet shape and divergence of 
the iteration procedure. The initial value of d is 0.5. During the itera- 
tion process, when the rate of change of the total free sheet vortex element 


force is less than five percent, d is set to 0.95. 
s the length of each segment 
As - Ax [l + (f£) 2 + iff) 2 ] 


Since the length of each segment is conserved, therefore, 

Az.2 1 4 


Hence 


(Ax) 


As 

[‘ * W 1 + 


It follows that 


(Ay)* - (Ax)* 0 t£) 


(Az) 


Ax 

(Ax) (^) 


(56) 

(57) 

(58) 

(59) 
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Thus, ths new location of ths (i41) th and point la glvan by 

x i+l " X 1 + 
y l+l " y l + (Ay) * 
z i+l " 2 i + (Az) * 


( 60 ) 

(61) 

(62) 


3.5.3 Formulation of Fores Frae Condition of Vortax Cora In Tha Cora Model 
During each iteration, the Joukowsky force acting on each segment of the 
core la computed. Let 
f - $ x ? 


- <ut + v} 4 w£) x (r t 4r J 4 r £) 

x y z 


Then 


r u 

X 


r x u 


r 

w _z _w ,Az> 

U • r x u " W 


_Z _ V 


v Ax ; u 


(64) 

(65) 


If the new position of the core is such that F ■ 0 and F_ ■ 0, then 

y z 

(-)' - (—)' - 0 (66) 
u Ax 


<£>’ - <?' ■ ° 


(67) 


Assume that 


(-) - (-) 
u u 

£) - (*)' 

u u 


Combining equations (64), (65) and (66), (67) yields the following equations 
for the i 6 * 1 segment : 


(£L)' . t*2L\ - (—*) 
l Ax'i u ; i 

X 


(M)' « (^2.) +( — 1) 

Mx ; i V Ax ; i u ; i 
,th 


The new slopes of the i segment of the concentrated core are taken to be: 

/&)* - (i _ x ) ( ^Z) + x (4^)' 

'‘Ax'! ' y M Ax ^ z v Ax'i 


(— )* 
'Ax' i 


< l ‘ x z H ll\ + x y Oi 


( 68 ) 

(69) 

(70) 

(71) 


Where X y and \ z are relaxation parameters and are discussed in Section 3.8. 


19 


Combining aquations (68) through (71) results In: 


" ( i _ Xy( T?\ 
( e>* “ + x 2 ( rJ ) 1 


(70) 

(71) 


During the Iteration process, when the rate of change of the vortex core force 
Is less than five percent, the relaxation parameters X y and X z of each core seg- 
ment will be further reduced by multiplying the factor 0.1. Thus, the move- 
ments of core segment will not cause excessive fluctuation in the aerodynamic 
characteristics. 


Since the segment length is conserved, it follows that 

,Az.2 


As - Ax [l + (f £) 2 + (||) 2 ] ! 
As 

[i ♦ (£>“ ♦ f 


(Ax)* - 


(74) 

(75) 


Hence , 


(Ay)* - (Ax)* (|£) 


(Az)* - (Ax)* (—)* 


Thus, the new position of (i+1) 
as follows 


th 


x i+1 - x x + (Ax) 

i 

y i+1 - y ± + (^) 

z i+i " z i + (Az)1 


(76) 

(77) 

end point of the concentrated core are found 

(78) 

(79) 

(80) 


In the present analysis, the effect of the secondary vortices has not been in- 
cluded. 


3.6 Centroid of the Leading-Edge Vortex Filament System 

After establishing the leading-edge vortex filament system based on exer 
cising Mehrotra's code (ref. 18) through two iterations, the centroid of that 
diffused vortex element system can be found as follows (ref. 20): 
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Consider a series of cross flow planes, procesdlng from the wing apex 
toward the trailing-edge. The centroids of the vortex filaments penetrating 


these planes are given by, 

- . ili r i y» 

J r ci 

- Ul r l 


J 


cj 


(81) 

(82) 


where, 


■ th 


, Zj are the centroid position in the j plane, 

f*H 

* z lj are t *' e l ntersecC i° n position of the 1 vortex filament with the 


jth plane, 


,th 


r t is the circulation around the 1 vortex filament, 

T c j is the concentrated line vortex strength at the j*"** plane, 
Oj is the number of vortex filaments penetrating the j** 1 plane. 


3.7 Modification of Z-Coordinates of the Initial Concentrated Core 

In order to have better starting solution so that the number of iterations 
and computer time can be reduced, the initial z-coordinates of the concentrated 
cqre from Section 3.6 are modified to take the experimental values summarized 
by Smith (ref. 4): 

a' < 1.2 
1.2 < a' < 2.2 


z ■ 0.154a' + 0.1 
v 


z - 0.1333a' + 0.158 
v 

z - 0.13a' + 0.132 
v 


(83) 


a' > 2.2 


where 


a' ■ Tan a /Tan 0 
0 - 90° - A 

z v is the nondimens ional z-coordinate of the concentrated core, referring to 
the local semi-span; A the leading-edge sweep angle; and a the angle of attack. 
Equation (83) will be used only in the simple wing configuration. 


The lnlelel configuration on the core model, on the right half of the 
wing, is shown in Flge. 6A through 6C. The top view of the initial free eheet 
and the modified y-coordlnate of the core location are shown on the x-y plane 
in Fig. 6A. The initial modified y-coordlnate of the core location is such that 
it is at least at a location equal to 60 percent of the local semi-span of the 
wing (see references 4 and 22). The rear view of the detailed free sheet shape 
and its relative vortex core location above the wing surface are shown on the 
y-z plane in Fig. 6B. The side view of the free sheet and the modified z-coor- 
dinate of the core location are shown on the x-z plane in Fig. 6C. Note that 
the initial modified z-coordlnate of the core location is based on equation (83). 


3.8 Formulation of Relaxation Parameters in the Core Model 

In the first two iterations of the core model, the core is forced to 
move either in the y or z direction. The direction depends on the sign and 
the difference between the initially computed lift coefficient and the reference 
lift coefficient based on the suction analogy (references 26 and 27). From 
the numerical experimentation, the relaxation parameters X y and X 2 , in equa- 
tions (70) through (73), are taken as follows: 

First iteration 

If the initially computed lift coefficient is smaller than the reference 
lift coefficient based on the suction analogy, the core will move down and 
the corresponding rel^gatlon^pajameter is 

X z 1 “ ” 1,5 (C L Tan a) (84) 

s 

If the initially computed lift coefficient is greater than the reference 
lift coefficient, the core will move out and the relaxation parameter is 

X 


X - -2.0 


(C. - C T ) Sin a 

L i L 


(85) 


where C_ is the reference lift coefficient, and is the initially computed 
s l 

lift coefficient. The initial set of X and X are designed so that the 

y l Z 1 
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movement of cor* segment* is proportional to th* difference in and . 

Also, the higher the angle of attack, th* larger the value of X or the lower 

y l 

the value of X is used. It was determined from numerical experimentation that 
2 1 

the initial core needs larger move-down movement in the z direction than the 
move-out movement in the y direction. 

Second iteration 

If the Initial movement of the concentrated core is in the z direction, the 

y-coordinates of the core will be adjusted next in the second iteration. The 

computed lift coefficient after the initial adjustment of the concentrated core 

is denoted by C^. The core will be caused to move out if is greater than 

and to move in if is less than . The corresponding relaxation para- 
s s 

meter is: 


(C T - C. ) Sin a 
L L s 


( 86 ) 


} 2 W L 

s 

On the other hand, if the initial movement of the concentrated core is in 

the y direction, the z-coordlnate of the core will be adjusted in the second 

Iteration. The core will move down if C. is less than C T and move up if C. 

L L s L 

is greater than . The corresponding relaxation parameter ls: 
s 


(c , - > 

V ■ °- 75 rd-ara - 
8 


(87) 


The relaxation parameters X and X have the same form as the initial set ex- 

y 2 z 2 

cept the magnitude is reduced by half. Special conditions are imposed for all 

configurations when the angle of attack is below 20 degrees. They are (1) 

the magnitude of X in equation (84) is to be reduced by half; and (2) for the 
Z 1 

angle of attack below 10 degrees, the angle of attack in equations (84) and (87) 

is set to 10 degrees to avoid the singularity in a in computing X and X 

Z 1 z 2 

Therefore, the movement of the core in the first two iterations for angles of 
attack below 20 degrees will not cause a large fluctuation in the aerodynamic 
characteristics. 
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The relaxation parameters X and X of each segment of the core are to be 

y * 

used from the third iteration on, based on the force gradients of core segments 
before and after the adjustment of the free sheet vortex elements. Thus, the 
relaxation parameters to the i** 1 segment of the core are: 

F - F 

VV'T^rt ’ <88) 

VV'Hr^K • IV" <8,) 

In the above, F^, F^ and F ^ are the segment forces of the core before 

and after the adjustment of the free sheet vortex elements, and IF is the 

algebraic sum of segment forces acting on the core after adjustment of free sheet 

vortex elements. The.bign of C and C ire determined so as to reduce the 

y i 2 i 

total force acting on the f ree sheet vortex elements. 

In the core model, the locations of free sheet vortex elements and concen- 
trated core have a very sensitive effect on wing loading. To avoid the div- 
ergence of the Iteration procedure: (1) the restrictions are Imposed on the 

movement of free sheet vortex elements and concentrated core as discussed in 
sections 3.5.2 and 3.5.3; (2) the end point af the first core segment near 
apex will remain fixed from the third iteration; and (3) the z-coordlnates of 
concentrated core are not allowed to be below z . given by equations (40 and (50), 
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3.9 Aerodynamic Characteristics 

In reference 26 , the more accurate expressions to calculate the overall 

aerodynamic characteristics are derived by Lan for the cambered wing. Both 
dz 

and Yy" are Included in the computation. 


3z 
c 

ax 


Assume that the wing surface (Fig. 7) can be described by: 


z - z c (x,y) 


(90) 


where z c (x,y) is the ordinate cf camber surface measured from the x-y plane. 


Introduce a function f such that: 


f - z - z c (x,y) 


(91) 


Therefore, a unit normal vector on the wing surface can be defined as: 


3z 3z 


(92) 


l$fl r 3z c 2 3z c 2 1 

11 L 1 + ( inr> + ^ J 

The unit vectors and t^ are related with n (Fig. 8) in the following expressions 

n ■ -Sin a t + Cos ot ft (93) 

00 

ft ■ Cos a t + Sin a ft (94) 

00 

The computed lifting pressure AC^ (see Appendix A) is acting normal to the local 
camber surface, l.e., in the n diiection. Then, this pressure force can be de- 
composed into n and ft directions to determine the lift and drag components, 


respectively. Hence, the component in the lift direction is: 
3z r 3z 2 ^z 2 1 4 

Sin “ + Cos a)/ L 1 + + ( 3y £) J 


(95) 


the component in the drag direction is: 

3z 

c 


3z 


3z 


n-t . ’ ( * 


r °‘c 2 "c 21 1 
Cos a + Sin a)/ |^1 + (-^p) + (-g~) J 

To find the pitching moment about y axis, the ft component of (Ax ft + Az ft) x n 
is needed, where Ax and Az are the moment arms of the pressure elements refer- 
enced to the coordinate system. In the present analysis, Az is assumed to be 
zero. Hence, the corresponding ft-component of (Ax ft x n) is: 




( 97 ) 
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3.9.1 Sectional Aerodynamic Characteristics 

To find the sectional characteristics, it Is assumed that the wing vor- 
tices are situated along the cambered surface. Since the resulting lifting 
pressure Is acting normal to the local camber surface, the sectional lift, 
Induced drag and pitching moment coefficients can be determined by combining 
the lifting pressure AC^ and equations (95) to (97) through the local chord 
Integration. Hence, the sectional lift, Induced drag and pitching moment co- 
efficients at the j*'*' strip have the following expression: 



*j V * 4 

j 

r 3z 

- ' tj t 


3z 

j AC ( - "t— ~ Cos a + Sin 

p dX 


Cos a + Sin 


a)/ [ 

a V £ 


3z c 2 3z c 2 V 

1 + + <ar> ] 


dx 


r 

3z„ f 3z 

■f i “>■«*>' [ 

1 + <5T> + <lf> ] 


3z . 3z 

1 + ( — — ) Z + ( — -) 
v 3x * i v 3y } 


dx 




(99) 


( 100 ) 


J 

where the c terms In equations (98) and (99) are the leading-edge suction 

J 

effect; x. and x the leading-and tralling-edges x-coordinates of the chord 

j 'j th 

passing through the control station of the j strip; c the mean geometric chord, 

Cj the local chord length. The subscript in equations (98) and (99), means 
3z 3z 

C C fV 

(g^-) and (■ 577 -), are evaluated at the leading-edge of the j ’ strip. Using 


j ‘ 3y J 


this transformation. 


x - x. + tA- (1 - Cos 9) 

J 


equation (98) can be reduced to: 


■f 

J o 


AC p Sin 9 (£ 


Sin a + Cos 


a)/[l + ( 


♦ <£, 2 f 


de 
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► 3z 

+ v fe 


) Sin a + Cos a 

j 


] / b* <£>j * 


NTJ. p OS 1 r OS , OS 

rarir £ 4C p k [‘nr’, S1 “ • * '»• °] S1 " vL 1 * ( >r\ * ( ir’ J 

[ 3e ■* p 3z a 3z a«i H 

Si " • + c °* •].. / 1 1 + <inr>* + 'nr’jjt. 


( 101 ) 


and 


9 . 12 k , - 1 > 
k 2(N + 1) 


k-1, (N + 1) 


where the mid-point trapezoidal rule has been used to reduce the integral in».o 
a finite sum (ref. 19). Similarly, the sectional induced drag and pitching 
moment coefficients for the J th strip are given by: 


«T1 p OZ -I m OZ 

WT-ryZ; ^ L- <ir\ c “ “ + S1 " “J S1 ” VL 1 * <nr> 


3z - 

( £\ 2 

'3y ; 

y j 


1 


( 102 ) 


nii 


2*(N 

3z . 

(— ) 2 

( 3y \ 


»v“A _ If. o* * 

X ^ [x + ^(1-Cos 0 )]Sin0^[l+(_£) 2 + 

N + 1) k=l P kl- *j 2 * J k 

r 


(103) 


According to equations (52) and (53) of reference 19, the sectional 
leading-edge thrust coefficient is given by: 

. 1(1 - M* Cos* A) 4 ["i “t 008 ° ~ Sln a) ] 2 


'j 


2N 2 Cos A (1 - M 2 + Tan 2 A) 


(104) 


where is the free stream Mach number; A the sweep angle of the leading-edge; 
dz 

wj and (j— ) are the induced normalwash and slope of the wing surface at the 


j 


j 


leading-edge. For sharp-edged wings, c fc is assumed zero. 
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3.9.2 Overall Aerodynamic Characteristic* 

The total Induced drag, lift, pitching moment and leading-edge thrust 
coefficients of the wing are detarnlned by spanvlse integration of the sectional 
characteristics. Again, the integration is first transformed to an angular 
coordinate (0 <_ <_ it), and then reduced to finite sum by the conventional 

trapezoidal rule (ref. 19). Therefore, the total lift coefficient has the 
following expression, 

* 8 fa' 2 C * C ^ <105> 


where b is the wing span and C the wing area. 

y ■ (1 - Cos t>) 

equation (105) can be reduced to: 

IT 

c Sin p d<p 



and 



M-l 

c # c Sin <i 
i-1 *1 1 


By the transformation: 


♦ t - fjir , i-1, (M-l) 


(106) 


(107) 


where M-l is the total number of spanwlae strips. 

Similarly, the induced drag, pitching moment and leading-edge thrust 
coefficients are given by: 


c rt a 

D i 

2SM c di ^ Sin *i 

(108) 

3° 

IK 

^?i Cm i Cl Sln #1 

(109) 

C at 

T “ 

m C i Sln *i 

(110) 


3.10 Summary of Solution Procedures 

In the diffused vortex filament model of Mehrotra (ref. 18), the basic 
unknowns of the problem are the bound vortex density on the wing, the strengths 
and locations of the elements of the leading-edge vortex system and the 
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c railing waka. Tha problaa ia aolvad in eha following learatlva manner: 

(a) Praacrlba tha vortex- let ties for eha wing aurfaca, and eha 
initial locatlona of eha leading-edge and aide-edge fraa 
element a ovar eha wing and in eha waka. 

(b) By eatlafylng eha wing boundary condition, aquation (37), ob- 
tain eha bound voreax danalty of eha wing and eha atrangtha of 
fraa elemanta. 

(c) Calculata eha aerodynamic charactarlatlca and loada. 

(d) Adjuat eha fraa alamanta of eha leading- and slda-adgaa voreax 
ayatam and the trailing waka in eha local valocley vactor dir- 
ection. 

(e) Repeat atepa (b) through (d) for two iterationa. 

After two iterationa in Mehrotra's diffuaad vortax filament model, the 
present model with the leading-edge vortex core la ehan introduced (Sea Fig. 9), 
In the core model, the basic unknowns of the problem are the bound vortex 
density on the wing, the strengths and locations of the free sheet vortex 
elements, the concentrated core and the trailing wake. The problem is still 
solved in the iterative manner: 

(f) Find the centroid of the established leading-edge vortex 
filament system, equations (81) and (82), from step (e). 

(g) Allow the leading-edge vortex elements to merge and feed 
vorticity into the concentrated core through connecting 
segments. 

(h) From the numerical experimentation, the z-coordinate of the 
initial core is modified by equation (83), and the y-coordinate 
of the initial core is adjusted so that the Initial core 
location is at least 60 percent of the local semi-span 

of the wing to have a better starting solution. 
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(i) By satisfying tha wing boundary condition, aquation (37), 
tha initial bound vortax danslty of tha wing and tha 
strangth of fraa vortax elements ara obtalnad. 

(j) In tha fl'/st two ltaratlons of cora nodal, tha cora is 

forcad to nova althar in tha y or z dlractlon. T*.,« dir- 

action dapands upon tha sign and tha diffaranca batwaan 

tha initially conputad lift coalflclant and tha rafartnca 

lift coafflclant basad on tha suction analogy. Tha relax- 

ation parameters A and A are basad on aquations (84) 

y * 

through (87). 

(k) Fron tha third itaratlon of tha cora modal, tha movamant of 
tha cora is basad on tha relaxation parameters for each seg- 
ment of the core, aquations (88) and (89). These are used 
to keep the core moving in a direction for which the total 
algebraic sum of the forces acting on the free sheet is to 
be decreased. 

(l) Based on step (k) , the orientation of the leading-edge free 
elements of the free vortex sheet, side-edge vortices, the 
concentrated core and the trailing wake are adjusted in the 
local velocity vector direction. 

(m) By satisfying the wing boundary condition, equation (37), 
w^e bound vortex density of the wing and the strength of 
free elements are obtained. 

(n) Calculate all the aerodynamic characteristics and loads. 

(o) Calculate tne force acting on all free elements of the free 
vortex sheet over the wing surface and the concentrated core. 

(p) Repeat steps (j) through (o) until a converged solution is 
obtained. 
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The converged positions end orlentatlone for the free sheet vortex 
elements t side-edge vortex elements, snd the concentrated core are assumed 
to have occurred when the sum of the magnitude of the total force acting on 
the core and free sheet has reached a local minimum value in the iteration 
cycle. 

Numerical results to date indicate that when the angle of attack la 
decreased below 20 degrees, the present method has difficulty in generating 
a satisfactory starting solution. The main reason is the poor initial shape 
of the free vortex sheet from Mehrotra's model. Thus, with the core model 
method, the starting solution for wings at angles of attack less than 20 de- 
grees is taken to be the initial shape for 20 degrees angle of attack. 
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4. Numerical Results and Dlacuaalons 

This aacclon prasanta numarlcal charactarlatlcs of and pradlctlona by cha 
praaant mathod. Thesa pradlctlona ara compared with raaulta from tha mathods 
of Mahrotra (raf. 18), Johnson, at al. (rafa. 13 and 14) (callad Fraa Vortex 
Shaat (FVS)), Lan (raf. 26), Lamar (raf. 27) (called Vortex-lattice Method baa 
ed on Suction Analogy (VLM-SA)) and with data. Numarlcal characteristics 
and results of the present method will be presented and dlscusaed next. 

4.1 Numerical Characteristics 

The numerical results predicted by the present method have the following 
numerical characteristics: 

(a) Convergence Characteristics 

In the present method, numerical convergence can be expected during 
the iteration process If certain restrictions are Imposed on the movements 
of the free sheet and vortex core (see sections 3.5 and 3.8), 

(b) Initial Vortex Shape 

For the angle of attack below 20 degrees, the starting solution is 
taken to be the initial vortex shape from the 20 degrees angle of attack 
solution. 

(c) Spanwlse Strips 

In reference 18, Mehrotra made a parametric study to find a relation 
between the aspect ratio and the number of spanwlse strips which should be 
used to represent the lifting pressure on delta wings in order to obtain 
reasonably accurate aerodynamic results. The relation just mentioned is re- 
produced and displayed graphically in Fig. 10. From this figure eight 
vortex strips over the wing semi-span are seen to be needed for simple 
deltas in the aspect ratio range from 1.0 to 1.75. Because of this study 
and because of the strake-wlng and double delta configurations to be treated 
here being also of low aspect ratio, the same number was used Initially. 
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Du* to th* l**dlng-*dg* kink th**« configuration* had to b* r*pr«**nt*d by 
two spanwia* region* par semi-span. Exp*rl*nc* with th* QVLM of Lan's 
(r*f. 19) has shown that for two spanwis* region solutions, it is batter 
to have a different number of spanwlse strips in each region. The final 
choice of spanwise strips used her* is four for the strake region and five 
for the wing region. 

(d) Chordwlse Vortex Elements 

All results to be presented have been calculated by using six chordwlse 
vortex elements for all configurations. 

(e) Reference Lift Coefficient 

The reference lift coefficient discussed previously (see section 3.8) 
for use in adjusting the core movement is based on the suction analogy as 
implemented by Lan (ref. 26). 

(f) Complete Leading-Edge Separation 

In this study, only wings having complete leading-edge separation will 
be considered. 

4.2 Numerical Results 

This section presents numerical results by the present method based 
on the preceding ideas. These results are compared with other methods 
and data for certain test configurations which are shown in Figs. 11 to 15. 
Note that they consist of a flat delta wing (Fig. 11) , a conical cambered 
delta wing (Fig. 12), a delta wing with leading-edge vortex flap (Fig. 13), 
a double delta wing (Fig. 14) and a strake wing (Fig. 15). The one vortex 
core system is used for the delta wings while the two vortex-core system is 
used for the others. 


Overall Forces and Moments 


The overall aerodynamic characteristics are presented in Figs. 16 to 20. 
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(a) Planar 

The results for the flat delta wing are shown In Fig. 16. It Is seen 

that C. and C n can be accurately predicted by the present method as well as $ 

* 

by the other methods. However, the present method predicts the to be i 

more positive (l.e. more nose up) at high angles of attack, whereas the i 

methods of VLM-SA and Lan predict more nose down pitching moment. This j 

nose up moment of the present method could be due to the estimated peak 

■4 

| 

load being too forward. Mehrotra's results have good agreement with data 

(ref. 22) even though he used a diffused vortex model. j 

(b) Conical camber j 

The results for a conical cambered delta wing are presented in Fig. 17. | 

The results from Lan's method have the best agreement with data. The 
results are underpredicted by the present method at lower angles of attack 
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and by Cha VLM-SA at hlghar angles of attack. FVS raaulta (raf. 30) ara 


ganerally low. Tha C-. raaulta of Fig. 17(b) ahow that the FVS method la 

beat at angles of attack below 20 dagraaa and the preeent method la best 

for angles of attack above 20 degrees. The other methods all predict 

values of C n which are too low over the angle of attack range. Similar 
D i 

to the case for the flat delta, the present method predicts a pitch up 
tendency for the cambered wing. 

(c) Leading-edge vortex flap 

Recently, there has been much interest in applying the leading-edge 
vortex flap to Increase the lift-drag ratio under maneuvering conditions. 

Some low speed results for a 74°-delta wing have been reported in reference 29. 
The geometry of the planform is shown in Fig. 13 and it was modeled theore- 
tically with its undeflected trailing-edge flap omitted for simplicity. The 
results are presented in Fig. 18 where the C^, and have been based 
on the total planform area excluding the leading-edge flap. It is seen 
that the predicted results by the methods of Lan and VLM-SA agree well wi;h 
data. The and results estimated by the presented method slightly 
exceed the values of the data at moderate angles of attack. The re- 
sults show the present method to have a more nose ito moment than the data 
for angles of attack above 12 degrees. 

Having compared the results in the one vortex-core system, it is of 
interest to see the aerodynamic characteristics in the two vortex-core system. 


Double Delta Wing 

For the double delta wing (ref. 28), the comparison is made in Fig. 19. 

The results show that the C. and C_ are well predicted by the present 

L 

method and the methods of Lan and VLM-SA. At high angles of attack the 
present method predicts the results better chan others. 
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Strake Wing 


From reference 28, it is understood that the aerodynamic advantage of 
atrake-ving planforme are (1) the transonic maneuvering capability, and 
(2) the utilization of vortex lift. For the atrake-ving configuration of 
Fig. 15 (ref. 31), the results are presented in Fig. 20. It is seen that 
the predicted by the present method has good agreement with data whereas 
the method of Lan overpredicts the data in the high angle of attack region. 
The experimental data in Fig. 20(b) shows tha,. the pitching moments become 
more positive as the angle of attack is increased. This could be due to the 
wing vortex breakdown on the main wing. From the geometry of strake wing in 
Fig. 15, it can be seen that the wing area is much larger than the strake. 

As shown in the data in Fig. 20(a), the values start to decrease at about 

14 degrees of angle of attack, signifying the start of the wing vortex 
bursting. 


Lifting Pressure Distribution 

Figs. 21 through 29 show the spanwlse pressure distributions for different 
planforms at specified angles of attack and for constant x- locations. The 
results are discussed below. 

Delta Wing 

(a) Planar 

For the flat delta wing, the spanwlse pressure distributions are cal- 
culated at angles of attack equal to 10.2°, 20.4° and 30.7°, and compared 
with data at the two x-locations, x/c R - 0.5333 and 0.9333 in Figs. 21 to 23. 

In general, the pressure peak obtained by using Meoretra's method is lower 
than the data near midchord, higher near the traillng-cdge and is shifted 
toward the root chord. One objective in the present method is to improve 
the pressure distribution of Mehrotra's method, so as to allow for the occur- 
rence of a sharper pressure peak. The predicted pressure distributions 
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of the present method more closely resemble the date than do those of Meh- 
rotra. At the low angle of attack (a ■ 10.2°), the present method predicts 
a pressure peak that is Inboard of the data at midchord and is lower near the 
tralllng-edge. This may be due to the free sheet being too large, since its 
initial shape and core location is taken from the 20-degree solution and the 
size of free sheet will not shrink in the Iteration process. Thus, the in- 
board pressure distribution are higher than those measured. 

(b) Conical camber 

The results of the estimated spanwlse pressure distribution for the con- 
ical cambered delta are presented in Figs. 24 through 26, and are compared 
with data (ref. 22) and the FVS solution (ref. 30). Again, the presented 
results are at the angles, of attack equal to 10.2°, 20.4° and 30.7° and 
comparisons are made at two x-locations, x/c R ■ 0,67 and 0.935. At low and 
moderate angles of attack, the present method predicts values of the pressure 
peak which are lower than data, whereas the FVS solution predicts higher peaks. 
At high angles of attack the spanwlse pressure distribution predicted by 
the present method has better agreement with data than does the FVS solution. 
There are two pressure peaks near the trailing-edge at all angles of attack. 

The Inboard peak is expected and is typical for flat delta configurations. 
However, the reason for the outboard peak is not conclusively known. Some 
possible reasons for its occurrence are as follows: (1) the wing vortex 
distribution is not located on the highly cambered surface; Instead, it is 
on a flat surface; (2) the large spanwlse camber slope term (3z/3y = -1.4402) 
near leading-edge is not Included in the boundary condition; (3) eight span- 
wise strips may not be sufficient to resolve the lifting pressures in the tip 
region of this conically cambered wing. 

(c) Leading-edge vortex flaps 

Spanwlse pressure distributions predicted by the present method at 24-de- 
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gr«« angle of attack and two x-locatlona, x/c R • 0.5197 and 0.8937, are pr«- 
aantad in Fig. 27 for a 74°-delta ving with laadlng-adga vortax flap. Since 
there are no data or other theoretical raeulta available, only the preaent 
theoretical reeulta are preeented on each figure. Fig. 27(a) ahowa that 

the preaaure peak la at about 60 percent of the local aemlapan and that 
matches the vortex core location. Near the tralllng-edge r.he Interaction 
of the vortax core and leading-edge vortex flap produces large lifting 
pressure In the vortex flap region as seen In Fig. 27(b). 

Double Delta Wing 

From Fig. 28, the epanwlse pressure distribution at 20-degree angle of 
attack at two x-locatlona, x/c R ■ 0.5686 and 0.8099 are presented. In 
Fig. 28(a), the results near middle root chord indicate that there are two 
pressure peaks present. The lower pressure peak is behind the highly 
swept inboard region with the higher pressure peak being on the outboard 
wing region. The two vortex cores induce large sidewash in the tip region 
which result in a large pressure peak there. Fig. 28(b) shows the lifting 
pressures near the tralllng-edge and from this figure it is seen that another 
pressure peak has been produced. This outboard peak may be due to the large 
sidewash Induced near the tip due to the interacting inboard-located vortex 
cores. The exact reasons for this require additional investigation. 

Strake Wing 

From Fig. 29, the spanwlse pressure distribution at an angle of attack of 

12 degrees at three constant chord stations, x/c ■ 0.4,0. 7 and 0.9, are 

presenter. In Fig. 29(a), the results near the kink region show that the lifting 

pressure are mainly from strake vortex. In Figs. 29(a) to 29(c), the predicted 
pressure distributions are generally higher than experimental values, the 

large pressure peak near the tip region of each figure is mainly due to the 
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effect of side-edge vortex ay at am. Near: tha trailing adga tha flow flald la 
domlna.ad by tha wing with moat of tha predicted lifting praaaura balng on 
tha wing. It ahould ba notad that at an angla of attack of about 14° , tha 
wing vortax braakdown may hava occurred (raf. 37). Tha effect of vortex 
breakdown la not modelled in the praaent method. 

For a alender delta wing, a typical computer CPU time by tha praaent 
method la 3230 aeconda in 7 lterationa with CDC Cyber 175 at Langley Reaearch 
Center. This is to be compared with 4540 aeconda in 10 lterationa by Mehrotra's 
model. For the double delta wing analyzed in this report, the CPU time is 
5620 seconds in 8 iterations as compared with 22,450 seconds in 8 Iterations 
by the Boeing code (Ref. 14). If the first two iterations in Mehrotra's model 
can be by-passed by directly assuming a reasonable initial configuration, the 
saving in CPU time is about 1200 seconds for simple delta configurations. 

4.3 Mean Square Error Computation. 

To check the accuracy of the results predicted by the present and other 
theoretical methods, the so-called mean square error (MSE) will be used. 

Mean square error is defined as: 

MSE - - d 1 ) 2 )/K ] 1/2 (111) 

where is the i th experimental data value, K the total number of data points, 

and d^ the Interpolated theoretical results at the location of data point d^. 

Thus, the typical results of mean square error of each theoretical method to the 

overall lift coefficient and lifting pressure distributions are as shown in the 

following: 


MSE for Overall Lift Coefficient 
(a) a flat delta wing (Fig. 16(a)) 

Present Theory 0.0209 
Lan's Method 0.0172 
VLM-SA 0.0045 
Mehrotra's Method 0.0242 


_r y y* -.1— e* - 


(b) a conically cambarad wing (Fig. 17(a)) 


Praaant Thaory 
Lan'a Mat hod 
VLM-SA 
FVS Mathod 

(c) a dalta wing with L. 
Praaant Theory 
Lan'a Method 
VLM-SA 


0.0407 

0.0237 

0.0423 

0.0470 

E. vortex flap (Fig. 18(a)) 
0.0405 
0.0378 
0.0628 


(d) a double delta wing (Fig. 19 (a)) 
Present Theory 0.0351 

Lan'a Method 0.0330 

VLM-SA 0.0534 


(e) a strake-wing configuration (Fig. 20 (a)) 
Present Theory 0.0724 

Lan's Method 0.1197 
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MSE for Lifting Pressure Distribution 


Configuration 


* /c R 

a, das. 

Present Method 

Mehrotra 

Flat 

Delta, 

Fig. 

2 1(a) 0.5333 

10.2 

0.1694 

0.1198 

II 

tt 

Fig. 

21(b) 

0.9333 

10.2 

0.1391 

0.1539 

M 

II 

Fig. 

22(a) 

0.5333 

20.4 

0.2642 

0.3804 

ff 

M 

Fig. 

22(b) 

0.9333 

20.4 

0.1060 

0.1906 

M 

ft 

Fig. 

23(a) 

0.5333 

30.7 

0.5011 

0.7773 

♦ 1 

ft 

Fig. 

23(b) 

0.9333 

30.7 

0.1496 

0.4111 

Configuration 


X/C R 

a , deg. 

Present Method 

FVS Method 

Delta with 
Camber , 

Conical 
Fig. 24(a) 

0.67 

10.2 

0.1540 

0.1744 


If 

Fig. 

24(b) 

0.935 

10,2 

0.1470 

0.2163 


If 

Fig. 

25(a) 

0.67 

20.4 

0.3146 

0.3027 


If 

Fig. 

25(b) 

0.935 

20.4 

0.1690 

0.1358 


If 

Fig. 

26(a) 

0.67 

30.7 

0.5090 

0.4868 


It 

Fig. 

26(b) 

0.935 

30.7 

0.1491 

0.3139 


For the flat delta wing, the mean square errors for the lifting 
pressure Indicate that the present theory is better than Mehrotra's model. 
For a delta wing with the conical camber, the present results are also good. 
However, the lower values of MSE of present results partially may be due 
to the appearance of an unexpected pressure peak near the tip region. 
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3. Concluding Remarks 

A fro* ohoot-vortox cor* nodal to predict pressure loadings end 
overall aerodynasilc characteristics of low ispect-ratlo wings with edge 
vortex separation has been developed, by Incorporating a vortex cor* In 
the free sheet of Mahrotra's method. In the present method, Mehrotra's 
model Is used In the first two iterations to establish the vortex fila- 
ment system. Based on Smith's empirical results, the Initial vortex 
core location Is adjusted to improve the starting solution. Introduc- 
tion of a vortex core in the flow model Increases the numerical diffi- 
culty to obtain a converged solution. The problem is solved by using 
different sets of relaxation parameter!* for and by Imposing certain 
restrictions on movements of the free sheet and the vortex core. For 
low to moderate angles of attack ( a < 20* ), the starting solution Is 
taken to be the initial shape of the 20-degree angle of attack solution. 

Comparison with available data indicates that: (1) the present 

method Is generally accurate In predicting the lift and induced drag coeff- 
icients but the predicted pitching moment Is too positive; (2) the 
spanwlse lifting pressure distributions estimated by the one vortex core 

solution of the present method are significantly better than the results 

of Mehrotra's method relative to the pressure peak values for the flat 

delta; (3) the two vortex-core system applied to the double delta and 

strake-wing produces overall aerodynamic characteristics which have good 

agreement with data except for the pitching moment; (4) the computer 

time for the present method is about two thirds of that of Mehrotra's 

method for a delta wing and one quarter of that of FVS method for a 

double delta wing configuration. 

In view of the improvements embodied in the present method and 
possible further refinement, it would be of interest to: 
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(1) model tha fuselage portion as a potantlal flow region, with 
tha aaaunptlon of attachad flow In tha fuaalaga raglon and a 
vary aaall thruat coafflclant balng assumed In that raglon; 

(2) axamlna tha fuaalaga affact on wing vortax flow; 

(3) axamlna tha canard affact on wing vortax flow; 

(4) compara tha chordwlaa praaaura dlatrlbutlon with available 
data to Improve tha pitching moment computation. 
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Appendix A 

An Improved Lifting Pressure Expression 


From reference 18, it is understood that the total lifting pressure, for 

the wing with leading-edge vortex element system, is contributed by wing 

bound elements, wing streamwlse vortices, leading-edge bound elements and 

the streamwlse vortices of the leading-edge vortex element system. Hence, 

C H t H 

the total lifting pressure at the l c bound element of the j strip control 
station has the following expression: 

AC - (AC )_ + (AC ) (A 

"j.i p j.i t p j.t B 

where (AC )- is the lifting pressure due to streamwlse vortices; (AC ) R 

P j , i P J , i 

the lifting pressure due to bound elements. 

A. 1 Derivation of (AC )_ 

p j.i t 

From Fig. Al, it Is clear that the lifting pressure due to streamwlse 

vortices of wing and leading-edge vortex element are calculated along the 

common edge y^ of two adjoined spanwise strips. First, the force acting on 

the chordwlse element of length Ax of the leading-edge vortex system is: 

F ± - p vl ( ? i v t ) Ax (A 

where p is the fluid density, the free stream velocity, the vortex 
tH 

strength of the j leading-edge vortex element, and v^ the sidewash at 

point i. Thus, the force acting at the i t ' 1 point on the common edge per 

unit dynamic pressure and length is: 

(-?-) - 2 T. v. (A 

v qAx'i j i 

where q is the dynamic pressure. Similarly, the force acting at the same 
point due to the right leg of the strip, per unit dynamic pressure and 


length, can be written as: 
F„ / 


‘qdbj ■ “ 2 v if 1 Tdx 



Figure Al.- The leading-edge bound vortex element and wing streamwise 
vortices at the common edge 
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where y *• the wing bound vortex density, x^ the leading-edge x-coordlnate of 


the trailing leg under conelderatlon. Uelng the transformation, 
x - x 4 + ^ (1 - Cos 0) 

the equation (a. 4) can be reduced to the form: 

F_ r 9 


(A. 5) 


( 


ySin 0 d0 


R_x f 1 

Ax\ " " V 1 C j J 0 

ffC, V. r 1-1 

a - -J— [X \ Sin 9 k + H Y, Sin 9 t J 


(A. 6) 


where c^ Is the local chord length of the common edge at , N the number of 
bound elements in chordwlse direction and 0^ ■ —^ ir . The Integral has 

Veen reduced to a finite sum through the conventional trapezoidal rule. 


For the left leg of the (j+l) t * 1 strip: 


F. rc . v . |- i-1 

’ N Lj£ \ Sln 9 k + 15 Ti S1 » «J 


Therefore, the force per unit dynamic pressure and per unit length at the i 
point along the common edge is: 

F, 


(A. 7) 
th 


(... T. ) « (——) + (_ —) + (— t_) 

'qAx* Vx* V qAx * Vx\ 


(A. 8) 


Thus, the lifting pressure at i C ^ point along the j 1 "* 1 common edge other than 
the last one has the following expression: 


< 4C p j l ’' I ' 2 W/ <V2 - *J> 


(A.9) 


When the lifting pressure is calculated at the last common edge in the span- 
wise direction, the equation has the following form: 

F„ 


(AC 


)' " fe) / - y.) 


J.i 


qAx 


J+l ’3 


(A. 10) 


where y ^j+1 an< * are t * ,e y“l° cat * on °f common edges. Note that 
equations (A. 9) and (A. 10) are evaluated at all endpoints of wing bound 
elements, i.e. along the common edge. Surface spline interpolation (ref. 32) 
is performed to obtain the lifting pressure at the control station of each 
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vortex atrip. Therefore, (AC ) , the lifting pressure due to etreamwiee 

P J.i T 

vortices at the l** 1 bound element of the j** 1 strip control station is obtained. 


A. 2 Derivation of (AC ) 

Pj ’ iB h 

From Fig. A2, the normal force per unit length acting at the i^ bound 

element of the j** 1 strip is given by: 



p ( Ui Y y - v t Y x ) ASj 
2 

p (u t y t Cos Yj Sin <1^) Ay^/Cos 

2 q (Uj ■ v ± Tan ij^) y ± Ay^ 


where , 

Yj. is the bound vortex density, 

and v^ are the x and y components of the induced velocities, 
is the sweep angle of the bound element, 

AVj is the width of the j*"* 1 strip, 

As^ is the length of bound element, 

Y and y are the streamwise and spanwise components of the bound vortex 
x y 

density. 

t h tli 

Thus, the lifting pressure due to the i bound element of the j strip 
is : 

V , 

(AC n ) " " 2(U 1 “ V i Tan V Y i (A,11) 

P j,i B q ay j 

For a cambered wing without dihedral, additional lifting pressure will 

be generated from the interaction of a free stream component with the stream- 

wise vortex density y (ref. 26), as shown in Fig. A3. Adding this component 

of (-2 y Sin a Sin ) to the lifting pressure due to the wing bound element, 
x y 

equation (A. 11) is modified to: 
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(A. 12) 


(AC« ) 
J.l B 


K - 


Cos 


+ Sin a Sin A ) Tsn tj> | Yj 


where ^ may be defined as the local dihedral angle. Tan ^(dz^/dy), and 
may be different everywhere on a cambered wing. The correction term In 
equation (A. 12) Is applied only to the wing bound element and leading-edge 
bound element so that the local dihedral effect are included. 

At positive a, normally the wing bound elements produce positive lift- 
ing pressure near the leading edge; while the leading-edge bound elements 
produce negative lifting pressure. Thus, the leading-edge Kutta condition 
can be satisfied when the net lifting pressure Is zero at- the leading- 
edge. In order to substract the lifting pressure induced by the leading- 
edge bound element from that produced by the wing bounded vortex elements, 
the following prodecures will be used: 

(1) Extrapolate lifting pressure due to the wing bound elements in 
each chordwise direction to obtain the lifting pressure at the 
location of the leading-edge bound element EF. (See Fig. 3) 

Using the Fourier series relation, the lifting pressure at any 
chordwise location is determined by the following equation: 

AC Sin 6 — a + a„ Cos 2.0 

p o > 1 


(A. 13) 


2-1 


where Sin0 is Included to eliminate the known square root singularity of the 
lifting pressure at the leading - and t railing-edges. 

Therefore, 


a - — t* AC Sin0 d0 
o n 1 p 

o r 

1 N 

■v ^ AC Sin 0, 
r ' N 2_ P. k 

k-1 K 


A 

a.- — AC Sin0 Cos20 d0 

X TT 0 p 

q _ (2k-l)tr 

°k 2N ’ ’ 


2 N 

TT '5^ AC Sin0, Cos29, 
N P.. k V 


(A. 14) 


k-1 
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( 2 ) 


where N la the number of wing bound elecente in the chordwlse 
direction. The Integrels for Fourier coefficients ere re- 
duced to finite sums through the mld-polnt trapezoidal rule. 

Subtract suction pressure induced by the leading-edge bound 
elements. First, the concentrated strength of the leading- 
edge vortex element is converted to the vortex density. Since, 

T ■ / y dx (A. 15) 

from the transformation, 
c. 

x ■ x^ + (1 - Cos 0) (A. 16) 

it can be obtained that 


Tj ■ ^ fy Sin 0 d0 (A. 17) 

It is assumed that the concentrated vorticity due to the leading- 

edge vortex system is distributed near the leading-edge only. Using 

the mid-point trapezoidal rule, equation (A. 17) becomes: 

it c. Sin 0. 

1 j 2(N + 1) Y j 


(A. 18) 


or, 

2(N + 1) r 

Y j " it c. Sin 0, (A. 19) 

J j 1 

where 0, - y/. Thus, the decrease of lifting pressure 
1 Z(N + l; 

at the leading-edge bound element of the j th strip is: 


(AC ). - -2 fu. - (v. Cos 4 + Sin a Sin 4) Tan ij/.] y. (A. 20) 

p decrease L j J y y j -l* e i 

where the subscript £e means that u ^ , v ^ and are evaluated 

at the leading-edge bound vortex element of the j th strip. Hence, 

the actual lifting pressure at the location of the leading-edge 

bound element of the j 1 "* 1 strip is: 
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(ac ) - [<ac ) + (ac ) 1 + <ar ) 

r l.i L ? l.i* p 3.i T J *• ?J < 


(A.21) 


'decrease 

The lifting pressuree defined by equation (A.l; at a discrete number of 
points are Fourier-analyzed through equation (A. 14) so that lifting pressures 
at any locations can be calculated. 
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Appendix B 

Correction to Sldewash Calculation In Dlecrete Vertex Method 


In developing the free vortex filament theory for calculating aerodynamic 
characterlatlca of wlnga with edge-separated vortex flow, the poeltion of the 
concentrated vortex core le adjuated In the direction of decreasing Kutta- 
Joukowsky force. The magnitude of the force, and hence the new poeltion of the 
vortex core, depends on the calculated Induced velocity components on the core. 
It is well known that the induced sldewash cannot be accurately calculated by 
the discrete vortex approximation at points close to the vortex sheet. There- 
fore, the vortex core location and the wing aerodynamic characteristics will 
not be accurately evaluated. A method to Improve the accuracy of calculating 
sldewash by discrete vortex approximation will be detailed in the following 
paragraphs. 


B. 1 A Test Case 

Assume a vortex distribution with constant density (I ■ u t) exist in a 
rectangular region, as shown in Fig. Bl. The Induced velocity vector at any 
point in space with position vector $, is then given by: 




2 (lL - 1) X u 

r i 


6 r 
47 f 


d s 


'8 


ii ) /2 f a M (C - z)T - (n - y)k 

^ 0 ° [(€ - x > 2 + 8 2 (n - y) 2 + 3 2 (? - z) 2 ] 372 


d£ dn (B.l) 


where 

6 2 - 1 - M 2 

R 2 - (C - x) 2 + 8 2 (n - y ) 2 + 8 2 (C - z) 2 
^ ■ (£ - x)t + (n - y)T + (c - z ) 2 

Since the vortex density u is assumed to be constant o\er the entire region, 
equation (B.l) can be Integrated (ref. 33) to give a closed-form expression 
for the sldewash V as: 
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'A discrete vortex 


n i n i+l 


x,5 


Figure Bl.- Vortex dietrlbution In a rectangular region 


> 

N t n 



\ / 

\ / 


\< 

^ Q(x,y,z) 


21 and 34 are vortex 
segments of a strip 


Figure B2.- Local panel geometry 
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V 



<| - y)(a - x) 

Ic -«![(.- X> 2 + 6 2 ^ - y) 2 + B 2 (c - z) 2 ]** 

V>A 5 x—n 

- s| [(• - x) 2 + 8 2 y 2 + B 2 (C * z ) 2 ] H 

(; - y) (x) 

- z | [x 2 + B 2 (| - y) 2 ♦ B* (c - z ) 2 ] 4 


- Tan" 1 „ <~*> __ 

U - x | [x 2 + B 2 y 2 + 6 2 (c - z) 2 ]** (B,2) 

Equation (B.2) will ba uaad to check tha accuracy of approxlmaca mathoda of cal- 
culating tha intagrala In aquation (B.l). 


B.2 Approxlmata Mathoda 


The chordwlae Integration In aquation (B.l) la performed exactly ao that 
the aldewaah component becomea: 


V - is 

- 

iL 

8 2 . b/2 


w(n) d n J 

a - x 



T 

V 

o 

(n - 

y) 2 + (c - z) 2 \ 

[(a - x) 2 + 8 2 (n 

-y) 2 





A, — 

X 

1 

(B.3) 

+ 

B 2 

(C 


1 

X 2 + B 2 (n - y) 2 

+ 8 2 (t - a) 2 ]* 5 / 


In conventional 

discretized approximation, equation (B. 

3) la reduced to: 




iL 

Z ■ _ 

\ 

“"k / 

a - x 


v - 

V ■ 

4ir 

k-l (n k 

- y) 

2 + « -z) 2 l [( 

X 

a - x) 2 + fi 2 (n k - 

l 

y) i 

+ 

e 2 

(c 

- z) 2 ] 4 + 

r 2 

L x 

+ 6 2 (n R - y) 2 + 

B 2 ft - z) 2 Jj 

(B.4) 


The results In equation (B.4) are equivalent to applying mid-point trapezoidal 
rule to evaluating the apanwise Integral In equation (B.3). However, when 
(C - z) la email, the Integrand In equation (B.3) contalna a second-order sing- 
ularity at n ■ y, so that equation (B.4) Is expected to be inaccurate. To 
improve equation (B.4) and yet keep the method as simple as possible, the 
following consideration Is applied. 
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For convanianca, let 



to keep the method simple so that it Is easily extendable to more general 


applications. A correction factor may be applied, as will be Indicated later 


(see equation (B.10)). 







Substituting aquation (B.8) Into aquation (2,7) results in: 


v M 


(C - z)6 2 f J \ £ (n k> 

»-k-l 


4ir 


- H 


- *S 


(n k - y) + <c - z) 
<*> (y) f Cy) An 1+1 


(n 


i+l 


- yV 


^Tan” 1 


' 1+1 


- y 


C ~ 8 


J + aJ^Uzi 

+ (C - Z> 2 i TT^T 

-i "i - H 


w(y) f (y) An i 
(n t - y> 2 + (c 

a 2 . 


- zY 


- Tan 


C- z | 


(B.9) 


The first ten in equation (B.9) represents exactly what would be given by 
the conventional discretized approximation (see equation (B.4)). Therefore, 
the remaining terms represent a correction C: 


c . 0 ( <c - .)6 2 ,/“ <y) £W L \ 

C l (n, - y?T 


w(y) f(y) An 1+1 


(c - z) 2 (n i+1 - y) 2 


>n l+l I 

+ <C - z) 2 J 

+ ^ {xan- ^ - Xan'^J <B. 

where G is a correction factor mentioned above and is taken to be 1.1 in the 
following numerical study. 


B. 3 General Applications 

To apply equation (B.10) to general situations involving such as leading- 
edge vortex sheet, trailing-vortex mutual interaction, etc., all geometric 
quantities involved in equation (B.10) should be evaluated based on local 
panel coordinate system. To illustrate this, consider two vortex segments 
2T and 34 in Fig. B2 and "P" is a point on a plane normal to two diagonal 
vectors and fj) and is obtained by projecting the control point Q along 
the normal vector. The vectors ^ and Tj can be calculated as: 

$1 - <x 3 - Xj^) t + (y 3 - y^) 1 + (« 3 - zp t (B.ll) 

$ 2 " ( x 4 - x 2^ ^ 4 ( y 4 “ y 2 ) T + < z 4 “ z 2^ ^ (B.12) 
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It follows that a normal vector N may be defined as: 

t * 

$ - x $ 2 - x 3 “ x i y 3 " y l 


x 3 " X 1 


x 4 - x 2 


z 3 - Z 1 


z 4 “ z 2 I 


- t [<y 3 - y 1 )(s 4 - z 2 ) - (z 3 - z 1 )(y 4 - y 2 ^j 

+ 1 [(z 3 - Zj^Xx^ - x 2 ) - (x 3 - XjXz^ - z 2 )j 

+ t [(x 3 - x 1 )(y 4 - y 2 ) - (y 3 - y 1 X* 4 - « 2 )J 

- N t + N v 1 + N, £ (B. 

A unit vector associated with $ can be defined as: 

"n ■ (B. 

m 

To determine the coordinate of "P", the distance PQ is needed. However, 
PQ is just the projection of vector 2Q onto the normal vector n : 

PQ = 2? • n 

- (x - x 2 ) n x + (y - y 2 ) n y + (z - z 2 ) n z (B. 

It follows that the coordinates of point P are: 
x p - (x - x 2 )n x 

y p - (y - y 2 )n y 

z - (z - z~)n„ (B. 

p 2' z 

To find the distances "x" and "a - x" in equation (B.10) (through f(y)), 

the unit vectors u^ and u 34 are needed: 

(x x - x 2 ) 1 + (y x - y 2 ) 1 + (z x - z 2 ) t 

“ 21 [(x x - x 2 ) 2 + (y t - y 2 ) 2 + (z x - z 2 ) 2 ]^ 

.* (x 4 - x 3 ) 1 + (y 4 - y 3 ) t + (z 4 - z 3 ) k 

34 [< x 4 - x 3 ) 2 + (y 4 - y 3 ) 2 + (z 4 - z-j) 2 ]* 1 

Then, the "x" distance in equation (B.10) can be replaced by: 

X R * <V “21 + V “34 5 (B 


(B.13) 


(B.14) 


(B.15) 


(B.16) 


(B.17) 


(B. 18) 


(B. 19) 






where 


*2p " <x p ’ x 2 ) * + (y p " y 2> ^ + <z p * z 2 ) * 


(B.20) 


$ 3p ■ (x p - x 3 ) t + (y p - y 3 ) J + (z p - z 3 ) £ (B.21) 

On the other hand, the "a - x" distance in equation (B.10) can be replaced by: 

*R * 11 <V “21 + V “34> < B ' 22 > 

where 

$ lp ■ (x x - x p ) t + <y x - y p ) 1 + (z x - z p ) t (B.23) 

$ 4p ’ (x 4 - x p > * + ^4 " y p > T + < z 4 “ z p > % (B.24) 

To determine "y - n^" and "n x+x - y" in equation (B.10), cross products 

•f — ► 

of S^p with U 2 X and Sj p with u^ are needed. Hence, "y - n^" will be replaced 
by: 

^R* l ? 2p x “2ll <B>25) 

Similarly, "n x+x - y" will be replaced by: 

y R"^3p x “34l (B ' 26) 

Finally, "? - z" in equation (B.10) must be determined. It is simply 

replaced by -PQ in equation (B.16): 

C - v - - [(x p - x ) 2 + (y p - y ) 2 + (z p - z ) 2 (B. 27) 

B.4 Numerical Results 

The locations of control points chosen for illustration are Indicated 
in Fig. B3. The results are compared in Figs. B4 through B7 for M»0. Except 
for the case where the control point is directly under one of the discrete 
vortices (i.e. p^int B) , all results with the derived correction terms "C" 
in equation (B.10) appear to be good. In applications to free vortoi.:-?-* lament 
program, the control points on the concentrated vortex core are chosen to be 
nearly at the middle of vortex strips formed by the leading-edge free vortex 
filaments. 
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Appendix C 

Determining Initial Core Location by Vector Analysis 

After exercising Mehrotra's code (ref. 18) through two iterations and 
establishing the entire leading-edge vortex element system, the centroid of 
the leading-edge vortex filament system is determined from equations (81) and 
(82). Therefore, the initial configuration of leading-edge vortex filament 
system and its centroid are as shown in Fig. Cl. 

Based on the vector analysis of the cross product, the procedures of 
finding the intersection between the leading-edge vortex filament system and 
its line of centroid are discussed in the following. 

Note that each leading-edge vortex filament consists of a series of 
short straight segments. It is desired to find the intersection point betweeu 
the leading-edge vortex element system and its line of centroid on the X-Y plane. 
Thus, only x and y coordinates of each segment will be considered in the cross 
product analysis. Assume that S 12 is a segment along the leading-edge vortex 
filament system and !> 34 is a segment along the line of centroid. Then, $^ 2 
and ? 34 have an intersection if the following conditions are satisfied: 

(1) Resulting vectors of S 31 x S 34 and S 32 x ? 34 are in the opposite 
direction. 

(2) Resulting vectors of S^ 3 x 1^ 2 and S>^ 4 x 3^,, are in the opposite 
direction. 

For example, in Fig. C2, by using t!ia right hand rule and cross product analysis, 
? 12 and S 34 have an intersection point. Similarly, $ 12 and $ 34 have no inter- 
sections in Fig. C3. 
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Figure Cl.- Leading-edge vortex filament system and its line of 
centroid 










Appendix D 

Correction* to Induced Velocltie* Due to Vortex Core 


It 1* clear from Appendix A that the lifting presaure at the wing bound 
element will not be predicted accurately, If. the computation of aidewaah due 
to the vortex core la not treated properly. In fact, a very large aidewaah 
would reault whenever the vortex core la too cloae to the wing aurface. To 
improve the aidewaah computation due to the vortex core, and hence, the 
lifting pressure and the aerodynamic characteristics, the computed sidewash 
due to the vortex core will be corrected as follows. 

In reference 34 'on page 592), Sir Horace Lamb showed that the circum- 
ferential velocity around a circle with radius r due to a line vortex of 
strength r is: 


v 


6 



-r 2 /4vtj 


(D.l) 


Later, Kutler in reference 35 defined the approximate radius as: 


r - 4vt_ 
c o 


Thus, the tangential velocity in equation (D.l) becomes: 


v 


0 


r 

2rrr 


(1 - e 



2 

c) 


(D.2) 


(D.3) 


In equf' ion (D.l) through (D.3), r is the radial distance from the center of 
the concentrated line vortex, T the circulation strength, r c the approximate 
core radius, the kinematic viscosity, and t Q the Initial time. 

From equation (D.3), after expanding the exponential term, the tangential 
velocity becomes: 


v 


0 


r 

2irr 



2 

c) 




r 

2irr 



(D.4) 
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Th« above aquation la obtained from the aeaumptlon that the ratio r/r c la 

r 

leas than one. In the preeent anaiyaia, the tarn of la to be regarded aa 
the computed aidevaah without correction, and all terma In equation (D.4) are 
defined aa followa. 

In Fig. Dl, the vortex core aegment 3^ ig on th « right wing and $^g on 
the left wing, 3^ * leading-edge vector. Point 5 la the wing bound 

element location where the Induced velocltlea due to vortex core aegmente 
$24 *nd are to be calculated. For the right core aegment 3^, the 
character latic core radiuo r £ la defined aa 3/4 of the perpendicular diatance 
from point 4 to line 3^2*' 


r 

c 




’12 


14 


I*. 


12 ' 


(D.5) 


The numerical factor, 3/4, in equation (D.5) la an aaaumed empirical value. 

The radial distance r of the wing bound element at point 5 to the vortex 
core segment 3^ is calculated as: 

'Thus, the radial diatance r is the perpendicular distance from point 5 to 3^. 

Based on the assumption in equation (D.4), the computed induced side- 
wash v* is to be modified only when the ratio r/r c is less than one. The 
corrected induced sidewash is: 

v - v’ (-£) 2 , r/r < 1.0 (D. 7) 

r c — 

c 

that is the location at which the Induced velocity is to be computed is inside 
the vortex core region. 

Similarly, at the right wing bound element at point 5, the sidewash due 
to the left core segment 3 _„ is calculated, based on the following r and r: 

JO c 
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(D.8) 


r c * 4 


m 


’19 


18 


I* 


19' 


r ■ 


It 


’78 
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It 


■I • H 


f 1*14 

1 


t 


12 


1^12 1 


34 


64 


78' 


U 


34' 


(D.9) 


Again, tha correction of tha computed sidewaah will be applied only to those po- 
ints inside the vortex core region, and the modified sidewaah has tha same 
expression as that of equation (D.7). Note that the right wing vortex core 
segment 3^ and the left wing vortex core segment 3 7g are geometrically 
similar. Therefore, in the actual computation, r and r due to 3, Q arc 
based on the right wing vortex core segment 3^ with the bound element loca- 
tion switched from point 5 on the right wing to its image point 6 on the left 
wing. 
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*—■ control station 


vortex strip 



Figure 2.- Strake wing geometry without leading-edge and side-edge 


vortex system 


74 









Figure 5.- Vortex segment geometry 
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Figure 6.- Mathematical Model of the initial free sheet and vortex core 



Figure 6.- Continued 




Free vortex sheet 



Figure 6.- Concluded. 
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Figure 10.- Variation of number of spanwise strips with 
aspect ratio for simple delta wings 
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Figure 16.- Continued. 




Figure 17.- Longitudinal Aerodynamic Characteristics of a AR = 1.1493 Conical Cambered Delta; 





Figure 17.- Concluded 









Figure 18.- Concluded. 



Figure 19.- Longitudinal Aerodynamic Characteristica of a Double ai-g . 





Figure 20.- Longitudinal Aerodynamic Characteristics of a Strake Uing; 


Data (ref. 31) 
Present Theory 





Figure 21.- Lifting pressure distribution for AR - 1.147 Flat Delta Wing at 
10.2 degree angle of attack; M * 0. 




Figure 21.- Concluded 



Figure 22.- Lifting pressure distribution for. AR * 1.147 Flat Delta Wiag at 
20.4 degree angle of attack; M * 0. 




Figure 22.- Concluded 


Data (ref. 22) 



Figure 23.- Lifting pressure distribution for AR - 1.147 Flat Delta Wing at 
30.7 degree angle of attack; M * 0. 



Figure 23.- Concluded 


Data (ref. 22) 



Figure 24.- Lifting pressure distribution for AB - 1.1493 Conical Canbered 
Delta at 10.2 degree angle of attack; M • 0. 


Data (ref. 22) 
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Figure 24.- Concluded 
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Figure 25.- Lifting pressure distribution for AR ■ 1.1493 Conical Caabered 
Delta »t 20 A degr°° ang1“ atrark; 1* » 0. 








Figure 26.- Concluded 




Figure 27.- Lifting pressure distribution for AR = 1.147 Delta Wing with a Leading-Edge 
Vortex Flap at 24 degree angle of attack; 6 - 30* , M - 0.2 . 





Figure 27.- Concluded. 





Figure 28.- Lifting pressure distribution for Double Delta Wing at 20 degree 
angle of attack: M = 0. 



Figure 28.- Concluded. 



Figure 29.- Lifting pressure distribution for a Strake Wing at 12 degree 




Figure 2°.- Continued. 



